Explosive percolation with multiple giant components 
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We generalize the random graph evolution process of Bohman, Frieze, and Wormald [T. Bohman, 
A. Frieze, and N. C. Wormald, Random Struct. Algorithms, 25, 432 (2004)]. Potential edges, 
sampled uniformly at random from the complete graph, are considered one at a time and either 
added to the graph or rejected provided that the fraction of accepted edges is never smaller than 
a decreasing function asymptotically approaching the value a — 1/2. We show that multiple giant 
components appear simultaneously in a strongly discontinuous percolation transition and remain 
distinct. Furthermore, tuning the value of a determines the number of such components with smaller 
a leading to an increasingly delayed and more explosive transition. The location of the critical point 
and strongly discontinuous nature are not affected if only edges which span components are sampled. 

PACS numbers: 64.60.ah, 64.60.aq, 89.75.Hc, 02.50.Ey 



The percolation phase transition models the onset of 
large-scale connectivity in lattices or networks, in sys- 
tems ranging from porous media, to resistor networks, 
to epidemic spreading . Percolation was considered 
a robust second-order transition until a variant with a 
choice between edges was shown to result in a seemingly 
discontinuous transition Q. Subsequent studies have 
shown similar results for scale-free networks 0. H|, lat- 
tices 0, [l(| , local cluster aggregation models [ll[ , single- 
edge addition models [12], Il3|. and models which control 
only the largest component [l4j . It seems a fundamental 
requirement that in the sub-critical regime the evolution 
mechanism produces many clusters which are relatively 
large, though sublinear, in size [Til . H3 , Il5j . Most re- 
cently, the notions of "strongly" versus "weakly" discon- 
tinuous transitions have been introduced [lj|, with the 
model studied in Q showing weakly discontinuous char- 
acteristics, while an idealized deterministic "most explo- 
sive" model [HI, E3 is strongly discontinuous. Here we 
analyze and extend a related model by Bohman, Frieze 
and Wormald (BFW) [l| , which predates the more recent 
work, and show that surprisingly [TIj], multiple stable 
giant components can coexist and that the percolation 
transition is strongly discontinuous. 

The "most explosive" deterministic process [i~5l - lr7| be- 
gins with n isolated nodes, with n set to a power of two for 
convenience. Edges that connect pairs of isolated nodes 
are added sequentially, creating components of size k — 2, 
until no isolated nodes remain. The cutoff k is then dou- 
bled and edges leading to components of size k = 4 are 
added sequentially, until all components have size k = 4. 
k is then doubled yet again and the process iterated. By 
the end of the phase k = n/2, only two components re- 
main, each with size n/2. The addition of the next edge 
connects those two components during which the size of 
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the largest component jumps in value by n/2. Letting t 
denote the number of edges added to the graph, we define 
the critical t c as the single edge who's addition produces 
the largest jump in size of the largest component denoted 
AC mol ; here t c = (n - 1), with AC max = n/2. 

The BFW model begins with a collection of n isolated 
nodes (with n any integer) and also proceeds in phases 
starting with k = 2. Edges are sampled one-at-a-time, 
uniformly at random from the complete graph. If an 
edge would lead to formation of a component of size less 
than or equal to k it is accepted. Otherwise the edge is 
rejected provided that the fraction of accepted edges is 
greater than or equal to a function g( k) = l/2 + (2k)~ 1 / 2 . 
If the accepted fraction is not sufficiently large, the phase 
is augmented to k+ 1 repeatedly until either the edge can 
be accommodated or g(k) decreases sufficiently that the 
edge can be rejected. Explicit details are given below. 
Asymptotically, lim^oo g(k) — a with a = 1/2. 

BFW established rigorous results whereby setting 
<?(200) = 1/2, all components are no larger than k — 
200 nodes (i.e., no giant component exists) when to = 
0.96689n edges out of 2m sequentially sampled random 
edges have been added to graph. They further estab- 
lish that a giant component must exist by the time 
to = c*n out of 2m sampled edges have been added, 
with c* e [0.9792,0.9793]. Yet, they did not analyze 
the details of the percolation transition. We show that 
their model leads to the simultaneous emergence of two 
giant components (each of size greater than 40% of all 
the nodes), and show analytically the stability of the two 
giants throughout the subsequent graph evolution. We 
then generalize the BFW model by allowing the asymp- 
totic fraction of accepted edges, a, to be a parameter and 
show that a determines the number of stable giant com- 
ponents that emerge and that, in general, smaller values 
of a lead to a more delayed and more explosive transition. 

Stating the BFW algorithm in detail, let k denote the 
stage and n the number of nodes. Let u denote the total 
number of edges sampled, A the set of accepted edges 
(initialized to A = 0), and t = \ A\ the number of accepted 
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edges. At each step u, an edge e„ is sampled uniformly 
at random from the complete graph generated by the n 
nodes, and the following algorithm iterated: 

Set I = maximum size component in A U {e„} 
if (I < k){ 

A^Au{e u } 
u <- u + 1 } 
else if (t/u < g(k)){ k^k + 1 } 
else { u 4— u + 1 } 

Thus while t/u < g{k), k is augmented repeatedly until 
either k becomes large enough that edge e u is accepted 
or g(k) decreases sufficiently that edge e u can be rejected 
at which point step u ends. Note g(k) = 1 requires that 
all edges be accepted, equivalent to Erdos-Renyi [19(. 

We numerically implement the BFW model, and mea- 
sure the fraction of nodes in both the largest and second 
largest component, denoted C\ and C2, as a function of 
edge density t/n. As shown in Fig. fflja), two giants ap- 
pear at the same critical point and remain distinct. To 
establish that the BFW model shows a seemingly discon- 
tinuous transition we apply the numerical method intro- 
duced in [6j. Let t l {n) denote the last accepted edge for 
which Cifi < rO and t\(n) the first accepted edge with 
Cin > An, where n is system size and 7 and A are pa- 
rameters. A 4 (7, A) — t\(n) — t l (n) denotes the number 
of accepted edges required between these two points. As 
shown inset to Fig. QJa), A l (j,A)/n is sublinear in n 
and t l (n)/n and t\(n)/n converge to same limiting value 
t c /n = 0.976 (Fig. [Tib)) for both C x and C 2 . 

The discontinuous nature is made more explicit in 
Fig. [2fa) showing AC max , the largest increase of the 
largest component due to addition of a single edge, ver- 
sus n (blue squares are BFW). AC max is independent of 
n (strongly discontinuous). Essentially the same value 
of AC m ax is always observed and results from the sec- 
ond and third components merging together to overtake 
what was previously the largest. The model studied in Q 
(PR) shows a decrease with n (weakly discontinuous), 
with scaling n -°- 065 as also recently observed in |13l Il6j. 

The key to coexisting multiple giant components is the 
high probability of sampling internal-cluster links in the 
super-critical region which, by definition, do not increase 
the component size. We formalize this by first intro- 
ducing a function P(k, t, n) defined as the probability of 
sampling a random link which leads to a component of 
size no larger than stage k at step t for system size n: 



(1) 



d+Cj<k/r, 



where Ci denotes the fraction of nodes in component i. 
The first term on the right-hand side is the probability 
of randomly sampling internal-cluster links in all compo- 
nents. The second term is the probability of sampling 
spanning-cluster links which lead to a component of size 
no larger than k. This is valid for any configuration in 
phase k. We also consider S(k,n), the the probability of 
sampling random links which lead to components of size 
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FIG. 1: The BFW model, (a) Ci and C2 versus edge density, 
t/n, for n = 10 6 , showing the emergence of two stable giant 
components. Inset: A(l/2, 0.5)/n, A(l/2, 0.35)/n for Ci and 
Ci versus n. (b) The lower and upper boundaries of A/n for 
both Ci and Ci, where data points are averages over 200 to 
2000 realizations and dashed lines are best fits. 



no larger than k if all possible spanning-cluster links are 
added in stage k. Thus 



S(k,n) = J2C? 



(2) 



(Note the values of the CVs in Eq. [T] and [5] can differ 
from each other.) For any specific stage k, it is easy 
to show that P(k,t,n) > S(k,n) since, if t increases, 
P{k, t, n) can only decrease or stay the same. More ex- 
plicitly, if an internal-cluster link is added then P(k, t, n) 
is invariant, while if a spanning-cluster link is added be- 
tween components i and j then the first term increases 
by {d + Cjf - (Cf + Cf) = 2C i C j and second term 
decreases by at least 2CiCj. (Additional decreases result 
if there exist components / satisfying C'i + Cj < k/n, but 
with {d + Cj) + Ci > k/n). 

Focusing now on the critical region, let k* denote the 
value of k at t c . Numerical results for a variety of sys- 
tem sizes show that at t c when n > 10 6 , k* /n ~ 0.570, 
C\ ~ 0.570 and C2 ~ 0.405 with error bars of or- 
der O(10 -4 ) obtained over 30 to 300 realizations depen- 
dent on n. Thus the remaining components have total 
size density J2i>3 Ci — 0.025. We can establish a uni- 
form lower bound on S(k, n) for all k > k* using the 
simple intuition that under the normalization condition 
Ci = 1, S(k,n) = Cf is minimized when the num- 
ber of components are as numerous as possible and of 
similar size. Given that J2i>3 Ci < C\ — C2, the lower 
bound on S(k,n) is if all small components connect to 
C 2 . Then P(k, t, n) > S{k, n) > C\ + (C 2 + £ s :> 3 Ci) 2 = 
0.570 2 + (0.405 + 0.025) 2 - 0.510, so for any stage 
k > k*{n), we have 



P(k,t,n) > a = 1/2. 



(3) 



So for k > k* , the expected fraction of accepted links 
approaches a positive value strictly larger than a. 

Having established that in expectation P(k,t,n) > a, 
for k > k* , we need to explicitly consider what happens if 
an edge connecting the two giant components is sampled 
in this regime. Here (C\ + C 2 ) > k/n > k*/n and, by 
definition, t/u > g(k). Consider the case when edge e u +\ 
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FIG. 2: (a) AC max , the biggest single edge increase in Ci, is 
independent of system size n for BFW and for the restricted 
BFW process, indicating these transitions are strongly discon- 
tinuous, (b) Evolution of the distribution of n(s), the fraction 
of components of size s, for the BFW model. 



FIG. 3: (a) Number of stable giant components of size greater 
than O.ln versus a for n = 10 s . (b) When a = 0.3, three giant 
components emerge simultaneously. Inset shows convergence 
of the upper and lower boundaries of A/n for Ci, Ci, C3, anal- 
ogous to Fig. [Hb). 



connects C\ and C2. If tj (u+1) > g(k) the edge is simply 
rejected. But if tj (u+1) < g(k) either k needs to increase 
until the edge is accommodated, or (as we show next) a 
small increase in k quickly leads to t/(u + 1) > g(k). 
Setting t/u to the smallest value possible: 




Differentiating both sides in Eq. Q by fc we find that 

du 1 nt 

d(k/n) ~ 2^2(1/2 + y/T0)^k^' 

After the critical point we know that t ~ 0(n) and the 
stage k ~ C\n ~ 0(n). Thus from Eq. ([5]) it follows that 
~ £>(n 1/2 ) asn^co implying that an ©(n^ 1 / 2 ) 

increase in k/n results in t/(u + 1) > g(k), so the link 
which would lead to merging C\ and C2 is rejected and 
the two giant components are stable throughout the sub- 
sequent evolution. We verify this numerically. Letting 
k(n) denote the largest value of the stage ever attained 
for system size n, we find (k(n) — k*(n))/n ~ n~ 7 with 
7 = 0.46 ±0.03, and as n — > 00, k*(n)/n and k(n)/n con- 
verge to the same limiting value of approximately 0.570. 

The BFW model samples edges uniformly at random 
from the complete graph. If we restrict the process to 
sampling only edges that span distinct clusters, we ob- 
serve that two components with the same C% = 0.570 and 
C2 = 0.405 values coexist for several edge additions be- 
fore merging together. When they do merge the largest 
jump in Ci, equal to the size of C2, occurs. This is a 
strongly discontinuous transition as shown in Fig. [2ja) 
(the red diamonds) with jump size equal to 0.405. 

We now generalize the BFW model so that g[k) — 
a + (2fc) _1 / 2 (i.e., the asymptotic fraction of accepted 
links is now a parameter a). For the unrestricted pro- 
cess (sampling from the complete graph) we find that 
a controls the number of stable giant components. Let 
N(a,m) denote the number of stable giant components 
with size larger than m which appear at the critical 
point and remain throughout the subsequent evolution. 
Fig- ISI a ) shows N(a, O.ln) versus a for the unrestricted 



process, with system size 10 6 and each data point aver- 
aged over 100 independent realizations (showing no fluc- 
tuations). As a first decreases from a = 1, N(a,0.1n) 
increases, going from one giant component to two at 
a = 0.511 ± 0.003. Then, once a < 0.11, N (a, O.ln) 
decreases. (Using a less stringent criteria that considers 
all macroscopic components Cin > cn where c > a 
"giant", then N(a,cn) actually continues increasing.) 

The same reasoning that applied to the original BFW 
model can be used here to show the stability of the mul- 
tiple giants. Once k > k*, in expectation P(k,t,n) > a. 
Likewise, once k > k* , d ^j n ^ ~ C(n 1 / 2 ), so k/n increases 
very slowly and the process frequently samples new links 
and rejects links that merge any two giants. For exam- 
ple, if a = 0.3, N(a,0.1n) = 3 with C x = 0.414, C 2 = 
0.321, C 3 = 0.265, so P(k,t,n) > Cl+Cl + Cj - 0.345 > 
a = 0.3 when k > k*(n). See Fig.^b) for details. 

In Fig. |4ja) we show the typical evolution for the un- 
restricted BFW process for various values of a in the 
regime where only one giant component emerges. We 
measure the scaling window A(7, A), as discussed earlier, 
and find that smaller a leads to a more "explosive" tran- 
sition in that A is larger and the scaling window shrinks 
more quickly. Explicitly, for a = 0.7,0.8,0.9 (and setting 
7 = 1/2), we find respectively that A — 0.9,0.8,0.7 and 
t c /n fa 0.915,0.862,0.780. This delayed and more explo- 
sive nature with smaller a is intuitive in that the more 
links are rejected at each stage, the longer one stays in 
that stage, resulting in more components of size C^n ~ k. 

Figure HJb) shows the analogous behavior for the re- 
stricted BFW process (where only edges that span com- 
ponents are considered). The delayed and more explo- 
sive nature of the transition with decreasing a is also 
observed here. We also note that the location of t c is not 
affected. For instance, for a = 0.3, 0.5, 0.7, 0.8, 0.9 we 
find that t c /n « 0.998,0.976,0.915,0.862,0.780 for both 
the restricted and unrestricted processes. 

The behavior of the restricted process can also be ex- 
plained via Eq. (fTJ. Here because intra-cluster links are 
not allowed, the first term on the right-hand side van- 
ishes. If the stage stops at some ko < n, then the second 
term on right side of Eq. (TT]) decreases to which makes 
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FIG. 4: The results of varying a on (a) the unrestricted 
BFW process in the regime where only one giant component 
emerges, and (b) the restricted BFW process. Both (a) and 
(b) show increasing delay and larger AC msl for smaller a. 

P(fco, t,n) < a and stage keeps on growing until the two 
giants merge together. The restriction on sampled links 
does not change the nature of the transition. We find 
that the general BFW model is is strongly discontinuous 
for all a £ (0,0.97], regardless of whether link-sampling 
is restricted or unrestricted. 

In summary, we have analyzed the critical behavior 
of the BFW model and find that two stable giant com- 
ponents emerge at the same critical point in a strongly 
discontinuous transition (Fig. [2fa)). If we restrict the 
sampled links to only spanning-cluster edges, multiple gi- 
ants coexist for a few moments until they merge together 
in a larger discontinuous jump (Fig.[2](a)) and ultimately 
only one giant component emerges. We further general- 
ize BFW by making the asymptotic fraction of accepted 
links a parameter a, and find that number of stable giant 
components increases while a decreases. 

The existence of multiple stable macroscopic compo- 
nents is surprising and unanticipated (l8j . and has not 
been previously observed in stochastic percolation. A 
model of cluster aggregation where largest clusters are 



occasionally "frozen" and prevented from growing, does 
lead to multiple giant components 20] , however the freez- 
ing is imposed on the system. Simple algorithms that 
generate multiple giant components may create a new 
range of applications. In addition to providing insight 
and a potential mechanism for controlling gel sizes dur- 
ing polymerization [20], they may be useful for creat- 
ing communication networks consisting of multiple large 
components operating on different frequencies or for an- 
alyzing epidemic infections simultaneously arising in dis- 
tinct, independent groups. The unrestricted process al- 
lows multiple links between nodes and self-loops. For 
finite systems it may be important to understand what 
happens when only links not yet added to the graph are 
sampled. In the asymptotic size limit there should be no 
difference as 0(n) edges are added from 0(n 2 ) edges. 

The nature of the transition observed in Q was re- 
cently analyzed using cluster aggregation models with 
choice, where a set of candidate edges are simultaneously 
inspected at each step [2l[. The mechanism here, in con- 
trast, inspects one edge individually at each time. As 
shown in Fig.^a) the models in [2l[ (labeled MR (m=2)) 
and in [6] (labeled PR) show weakly discontinuous tran- 
sitions, where AC max decreases with system size. In con- 
trast both the restricted and unrestricted BFW models 
are strongly discontinuous, with a jump independent of 
system size. Finally we show the evolution of the com- 
ponent size distribution for the original BFW model in 
Fig. [U(b) (n(s) is the number of components of size s 
divided by n). This bi- modal distribution has a large 
right-hand tail, which deviates from a power law. 
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